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MAPPING CLOSURE APPROXIMATION TO CONDITIONAL DISSIPATION RATE 

FOR TURBULENT SCALAR MIXING 

GUOWEI HE* AND R. RUBINSTEIN 1 


Abstract. A novel mapping closure approximation (MCA) technique is developed to construct a model 
for the conditional dissipation rate (CDR) of a scalar in homogeneous turbulence. It is shown that the CDR 
model from amplitude mapping closure is incorrect in asymptotic behavior for unsymmetric binary mixing. 
The correct asymptotic behavior can be described by the CDR model formulated by the MCA technique. 
The MCA approach is outlined for constructing successive approximation to probability density function 
(PDF) and conditional moment. 

Key words, turbulent mixing, mapping closure, conditional dissipation rate, probability density function 
Subject classification. Fluid Mechanics 

1. Introduction. The probability density function (PDF) approach has been shown to be a useful 
tool in turbulence research. The systematic approach for determining the PDFs is by means of solving the 
transport equations for the PDFs. In the PDF equations for turbulent scalar fields, conditional dissipation 
rate (CDR) appears as the only unclosed term [1]. The recently developed large-eddy simulation schemes for 
turbulent reactive flow, such as filtered PDF approach [2], conditional moment closure [3] and Lagrangian 
fiamelet model [4] also appeal the CDR models. 

No satisfactory closure model for CDR had been constructed until the mapping closure approach [5, 6] 
was formulated. Amplitude mapping closure suggests a CDR model [7] whose form is separable in scalar 
and time variables. The model is in good agreement with direct numerical simulation (DNS) for initially 
symmetric binary mixing but fails in describing asymptotic behavior of the CDR for initially unsymmetric 
binary mixing [8]. Girimaji [9] has developed a novel amplitude mapping closure approach in which the 
reference fields are time-dependent. The CDR model obtained from a time-evolving Gaussian reference field 
still fails to describe the asymptotic behavior, but the one from a time-evolving Beta reference field can 
successfully describe the asymptotic behavior. This strongly suggests that the amplitude mapping closure is 
inadequate to describe the asymptotic behavior by itself. 

In this paper, we will develop a novel mapping closure approximation (MCA) to make successive ap- 
proximations to statistics of a scalar in homogeneous turbulence. This technique will be used to construct 
a CDR model which accounts for the asymptotic behavior of the CDR. In Section 2, we will investigate the 
asymptotic behavior of the CDR model from amplitude mapping closure and explain the reason why it fails 
to describe the asymptotic behavior correctly. In section 3, w r e will outline the MCA technique for successive 
approximation. In section 4, we will use the MCA technique to formulate a novel CDR model and compare 
it with DNS results. We will conclude with a summary in section 5. 
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2 . Asymptotic behaviors of amplitude mapping closure. In this section, we will show that the 
CDR model from amplitude mapping closure has incorrect asymptotic behavior for unsymmetric binary 
mixing. We consider the simple case of a single conserved scalar Z(x,t) in incompressible homogeneous 
and isotropic turbulence. With u(x,t) being the velocity field and k the molecular diffusivity, the evolution 
equation for scalar Z is 


(2.1) |^ + u VZ=kV 2 Z. 

Here the boundary condition is periodic and the scalar is initially homogeneous and isotropic. Therefore, 
the scalar Z(x, t) will remain homogeneous and isotropic for all time. We further assume that the scalar is 
binary (0 or 1) and its initial PDF is a double-delta distribution: 


( 2 . 2 ) 


P{Z,0) = { 


A 

l- A 

0 


if Z = 0 
if Z = 1 
if Z 0 or 1, 


where 0 < A < 1. .4 = 0.5 implies that the binary scalar has the same -weights and the corresponding PDF 
P{Z, 0) is symmetric with respect to Z = 0.5, while .4 ^ 0.5 implies that the binary scalar has different 
weights and the corresponding PDF P(Z, 0) is unsymmetric. DNS [10] has shown that the scalar PDF 
P{Z, t ) asymptotically approaches a Gaussian distribution, whether the initial double-delta distribution is 
symmetric or unsymmetric. The mean (Z) — 1 - .4 remains unchanged in turbulent mixing. The CDR 
X(Z, t) = k((VZ) 2 \Z) is nearly parabolic. At early stages, the CDR maximum is located at Z = 0.5, the. 
mean of the initial interface of the binary scalar. The maximum then moves and finally approaches the mean 
(Z) of the scalar, accompanied by a distortion of the parabola. Therefore, for an initially symmetric binary 
scalar (.4 = 0.5), the maximum will remain fixed at Z - 0.5, where the mean of the scalar and the mean of 
the initial interface are coincident. For an initially unsymmetric binary scalar (.4 ^ 0.5), the maximum will 
shift to the mean Z = 1 — A of the scalar from the mean Z = 0.5 of the initial interface. 

Amplitude mapping closure assumes 


(2.3) Z = X{<t> 0 ,t), 

where <&>(x) is a homogeneous Gaussian random field. The governing equation for the mapping function X 
is 


ox . 9 . ( <t > o ox o-x\ 

(2 ' 4) dt - ( <$> d<i>0 + a<& ) ’ 

where (j>o x is the spatial gradient of the scalar (f>o . The initial mapping function corresponding to ( 2 . 2 ) is 


(2.5) 

where 7 — y/2erf~ 1 {2A - 1 ). 


•V(4,0)= “ 

[1 if 00 > 7 , 

The exact solution of the equation (2 


.4) with the initial condition (2.5) is 
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( 2 - 6 ) 


x(<fo,t) 




where 


(2.7) 


^0 “ 


<fo 

v/1’ 


T — K 


(O 

(<AE> 




E = \/exp(2r) — 1. 


We can calculate the CDR from the exact solution (2.6) 


* (Z,() ' K wvhz exv ( _2 [frr ' (2Z ~ ‘fl 2 ) 

(2.8) = \(0.5, t) exp (-2 [erf~ l (2Z - l)] 2 ) . 

The most striking feature of the CDR model (2.8) is the separability of its form in Z and t. The 
separability suggests that the CDR's shape remains unchanged although its amplitude decays with time. The 
CDR maximum is located permanently at Z — 0.5 without any shift for either symmetric or unsymmetric 
initial conditions. The permanent location of the maximum for unsymmetric initial condition [7, 8, 9] is 
contrary to the known asymptotic behavior of the CDR. Therefore, the CDR model (2.8) from amplitude 
mapping closure is not able to describe the asymptotic behavior of the CDR. The amplitude mapping closure 
with a time-evolving Gaussian reference field gives a similar CDR model whose form is separable in t and Z 
and also fails to describe the shift of the maximum. Therefore, the separation in form is the main reason for 
the incorrect asymptotic behavior of the CDR model from amplitude mapping closure. 

3. Mapping closure approximation. The basic idea of mapping closure is to represent an unknown 
random field by mapping of a known random reference field. The statistics of the unknown random field 
can be calculated from the mapping function and the known reference field. The governing equation for 
the mapping function can be obtained from the evolution equation of the unknown random field and the 
governing equation of its PDF. 

Amplitude mapping closure assumes that the unknown random field can be mapped by a single known 
Gaussian reference field: Z = and a coordinate transformation dx'/dx = J(d<p(x , )/dx\t) that 

accounts for turbulent stretching. The assumption [11] holds if and only if the spatial level crossing frequency 
at which the unknown random field passes through a given value has a single maximum as a function of 
that value. Therefore, the existence of the mapping function as well as the coordinate transformation is not 
ensured for arbitrary unknowm random fields. Physically, a turbulent field exhibits eddies of different length 
and time scales so that it cannot be mapped by a single reference Gaussian field with a compact spectrum. 
Therefore, it is necessary to introduce more reference fields accounting for different eddies of different time 
and length scales: 


(3.1) 


Z = X{<f>o(x,t),<pi{x,t),--- t);x,t) 
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where 0*, i = 1 , — , ti , are reference fields. So far, we have not imposed any constraints on the reference 
fields so that we have some freedom in choosing the reference fields. For example, 0 O could be Gaussian but 
<f>\ could be a Beta random field. 

A one-to-one mapping could be established artificially between (Z, 0i, • * • ,4> n ) and (0 O , 4 >\ , * ■ * 

Thus, we can calculate the PDF and the conditional dissipation rate of the scalar 

(3.2) P(Z;xJ) = J (J^j d & 

(3.3) ((VZ) 2 \Z) = {(VX) 2 \Z). 

Here the integration is taken over the entire subspace of the composition 4>i The ensemble average 
is taken over the level surface on which (<£ 0 , <f >\ , • * , <t> n ) satisfies the constraint for a given Z 


(3.4) = Z . 

Introduction of more reference fields is expected to improve the approximation accuracy of mapping 
closure. The reference fields might be chosen to be statistically orthogonal so that more information can be 
introduced at less expenses. The hierarchy of mapping closure approximations with an increasing number 
of reference fields remains essentially non-perturbative. A perturbative approach, like the Wiener- Hermite 
expansion [13], could lead to the same result, but would require a prohibitively large number of terms to 
achieve the accuracy attainable with mapping closure based on a small number of reference fields. 

4. The model for conditional dissipation rate. Broad classes of mapping function are admissible 
to MCA approach. The form to be considered here is 


(4.1) Z = K(<fo(x),^i(*);0> 

where 4>q(x) is a homogeneous random Gaussian field in space and <£i (t) an inhomogeneous random Gaussian 
field in time. Amplitude mapping closure requires (t) = 0, which fails in the asymptotic behavior of 
conditional dissipation rate due to lack of an independent time-evolving reference field. 

Following the standard method [11, 12], we can formulate the mapping equation: 

(4.2) f ■ + ^(^iz> = -<»vn Z > +K ( v >m- 

The conditional averages in (4.2) can be evaluated by homogeneity of the velocity and scalar fields and 
Gaussianity of the reference fields 


(4.3) 


(uVY\Z) = 0, 


(M 7) - 
( dt ]Z) 2 ( 0 ?) 

(V 2 F|Z) = -00 


Ml 

dt ’ 

( 08 *) dr 

(0 o) d0o 


+ <08* > 


d 2 F 

00o 
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Fig. 4.1. Normalized CDRs \(Z, t)/F(t) for the initially unsymmetric double-delta distribution with A = O.G. Line (DNS) 
and plus (model): initial CDRs; Dash (DNS) and diamond (model): final CDRs . 


An exact solution can be obtained from (4.2) with the evaluated conditional averages (4.3) and the initial 
condition (2.5) which requires <f > i (0) = 0 


(4.4) 


Y = 


1 + erf 


<Po - 7<? r 

v/2S 


+ 1 - exp[^>i - ($i) + 




(2(Z)-1) 


We will calculate the CDR \(Z, t ) = /f((\T) 2 jZ), where ensemble average is taken over the level surface 
on which i) satisfies Y(<f>o, <f>l '• i) = Z. The expression obtained is then simplified by the first-order 

approximation of its Taylor expansion with respect to the mean (<p\ (t)). This leads to 


(4.5) 

where 


X(Z, t ) = F(f) exp 


-2 erf~ l [(2Z- 1) — (1 — 


\/ 1 + /</< i f 


mz) - 1 )] 


(4.6) 


F(t) = 


Uo ^ 


/« 


Vo 


v/2 iS’ 

<$> ’ 


/♦, = 


The time scales f<p 0 and are parameters of the present closure which must be input externally. Usually, 
they are reset at each time using a dynamical scheme, such as f 9o = k((VZ) 2 )/(Z 2 ) and = /^ 0 , so as to 
give correct evolution of the CDR. The CDR model (4.5) is compared with DNS of the diffusion equation 
with k = 1. In the case of homogeneity, the use of the diffusion equation to validate the CDR model for 
turbulent mixing has a reasonable justification. Figure 4.1 shows the initial and final CDRs obtained by 
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the DNS and the CDR model (4.5) for the initially unsymmetric binary scalar with mean {Z) = 0.4. For 
exhibiting the shift of the maximum, the CDRs are normalized by their amplitudes F(t) and (4.5) is rescaled 
by the scalar’s variances. 

The CDR model (4.5) is no longer separable in Z and t. Its shape will shift while its amplitude decays 
with time. It is easy to verify x(0.5, 0)/F(f) = 1 and *((Z),oc )/F(t) = 1 so that the location of maximum 
dissipation shifts from initial interface Z = 0.5 to the scalar’s mean Z = (Z). While (Z) = 0.5, the equation 
(4.5) recovers the CDR model (2.8) from the amplitude mapping closure. Therefore, the CDR model (4.5) 
correctly describes asymptotic behaviors of CDR: the CDR’s maximum asymptotically approaches Z = (Z), 
while initially being at Z — 0.5. We emphasize that (4.5) is the first term of the Taylor expansion of the 
CDR obtained from (4.4) with respect to the mean This simplified form gives the correct initial 

and asymptotic behavior of the unsymmetric binary mixing at the expense of an inaccurate description of 
the intermediate behavior. Indeed, it is well-known [9] that at short times, the CDR remains symmetric 
about the location of the initial interface. The symmetry point (location of maximum dissipation) begins to 
shift only after the elapse of a finite amount of time, depending on the extent of the initial asymmetry [14]. 
The present model predicts that the symmetry point begins to move immediately. A model which is also 
accurate at short times can be formulated from (4.4) in the general form 


(4.7) 


X(Z, t ) = F(t) exp 


-2 (or/" 1 [(2Z - 1) - (1 - f(t))(2(Z) - l)]) 2 ] , 


where f(t) is a time-dependent function with /(0) = 1 and /(+oo) = 0. The time- dependent function f(t) 
is responsible for the shift of the CDR’s shape. However, it is very difficult to write its expression explicitly, 
which needs a time scale as same as the case of amplitude function F(t). 

It has been shown that the mapping (4.1) is an appropriate approximation for a scalar gradient field. 
However, it is not expected that the mapping (4.1) makes the same accurate approximation to the scalar field 
as it does to the scalar gradient field itself. The reason is that a homogeneous random field is statistically 
orthogonal to its gradient field. The amplitude mapping (2.3) is an appropriate approximation to scalar 
fields but fails to approximate its gradient field. Generally speaking, the mapping closure carried out at the 
level of single-point PDFs is not valid for two-point PDFs such as gradient fields. This is the motivation to 
develop mapping closure approximation of higher orders. We expect that a higher-order approximation to a 
random field could produce a correct description of its gradient field. 


5. Conclusion. We have developed a novel model of the conditional dissipation rate for turbulent mix- 
ing. The model is able to describe the asymptotic behavior of the CDR for either symmetric or unsymmetric 
initial double-delta distributions. The amplitude mapping closure has unsatisfactory asymptotic properties 
for the shift of the CDR’s maximum. The problem can be solved using the extended mapping closure ap- 
proximation developed in this paper. Further research is to account for the evolution of scalar variance by 
extending the CDR model (4.5) . MCA can incorporate the effect of multiple time and length scales of 
practical interest in its predictions using more than one reference field. It provides a useful approach to 
describe statistics of turbulent mixing. 

The author (G. W. He) is grateful for the very valuable discussions on the conditional dissipation rate 
with G. Kosaly and H. Pitsch. Special thanks are due to G. Kosaly for pointing out the exact solution (2.6). 
We thank S. S. Girimaji for his valuable comments on the shift process of the symmetric point. 
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